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A memory usage scheme 1or performing wavelet processing 



(57) A system comprising a memory and wavelet 
processing logic is described. The memory is sized to 
include lines to store a band of an image and additional 
lines. The wavelet processing logic comprises a wavelet 
transform and access logic. The wavelet transform gen- 
erates coefficients when applied to data in the memory. 
The access logic reads data from the memory into the 



line buffers to supply data stored in the memory to the 
wavelet transform and to store coefficients in the mem- 
ory, such that after data stored at a first pair of lines is 
read from memory into the buffers of the access logic. 
The access logic reuses the first pair of lines to store 
coefficients generated by the wavelet transform that are 
associated with a second pair of lines different from the 
first pair of lines. 
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Description 

FIELD OF THE INVENTION 

5 [0001] The present invention relates to the field of compression and decompression; more particularly, the present 
invention relates to memory usage for performing wavelet processing. 

BACKGROUND OF THE INVENTION 

w [0002] The new JPEG 2000 decoding standard (ITU-T Rec.T.800/ISO/IEC 1 54441 :2000 JPEG 2000 Image Coding 
System) provides a new coding scheme and codestream definition for images. Although the JPEG 2000 standard is 
a decoding standard, the JPEG 2000 specifies encoding and decoding by defining what a decoder must do. Under the 
JPEG 2000 Standard, each image is divided into one or more rectangular tiles. If there is more than one tile, the tiling 
of the image creates tile-components that can be extracted or decoded independently of each other. Tile-components 

75 comprise all of the samples of a given component in a tile. An image may have multiple components. Each of such 
components comprises a two-dimensional array of samples. For example, a color image might have red, green and 
blue components. 

[0003] Attertiling of an image, the tile-components may be decomposed into different decomposition levels using a 
wavelet transformation. These decomposition levels conlain a number ol subbands populated with coeff.c.ents that 

20 describe the horizontal and vertical spatial frequency characteristics of the original tile-components. The coefficients 
provide frequency information about a local area, rather than across the entire image. That is, a small number of 
coefficients completely describe a single sample. A decomposition level is related to the next decomposition level by 
a spatial factor of two, such that each successive decomposition level of the subbands has approximately half the 
horizontal resolution and half the vertical resolution of the previous decomposition level. 

25 [00041 Although there are as manv coefficients as there are samples, the information content tends to be concentrated 
in just a few coefficients. Through quantization, the information content of a large number of coefficients is further 
reduced. Additional processing by an entropy coder reduces the number of bits required to represent these quantized 
coefficients, sometimes significantly compared to the original image. 

[0005] The individual subbands of a tile-component are further divided into code-blocks. These code blocks can be 
30 grouped into partitions. These rectangular arrays of coefficients can be extracted independently. T-he individual bit- 
planes of the coefficients in a code-block are entropy coded with three coding passes. Each of these coding passes 
collects contextual information about the bit-plane compressed image data. 

[0006] The bit stream compressed image data created from these coding passes is grouped in layers. Layers are 
arbitrary groupings of successive coding passes from code-blocks. Although there is great flexibility in layering, the 
35 premise is that each successive layer contributes to a higher quality image. Subband coefficients at each resolution 
level are partitioned into rectangular areas called precincts. 

[0007] Packets are a fundamental unit of the compressed codestream. A packet contains compressed image data 
from one layer of a precinct of one resolution level of one tile-component. These packets are placed in a defined order 
in the codestream. 

40 [0008] The codestream relating to a tile, organized in packets, are arranged in one, or more, tile-parts. A tile-part 
header comprised of a series of markers and marker segments, or tags, contains information about the various mech- 
anisms and coding styles that are needed to locate, extract, decode, and reconstruct every t.le-component. At the 
beginning of the entire codestream is a main header, comprised of markers and marker segments, that oflers similar 
information as well as information about the original image. 

is [0009] The codestream is optionally wrapped in a file format that allows applications to interpret the meaning of, and 
other information about, the image. The file format may contain data besides the codestream. 

[0010] The decoding of a JPEG 2000 codestream is performed by reversing the order of the encoding steps. Figure 
1 is a block diagram of the JPEG 2000 standard decoding scheme that operates on a compressed image data 
codestream Referring to Figure 1 , a bitstream initially is received by data ordering block 101 that regroups layers and 
so subband coefficients. Arithmetic coder 1 02 uses contextual infonnation collected during encoding about the bit-plane 
compressed image data, and its internal state, to decode a compressed bit stream. 

[0011] After arithmetic decoding, the coefficients undergo bit modeling in coefficient bit modeling block 103. Next, 
the codestream is quantized by quantization block 104, which may be quantizing based on a region of interest (ROI) 
a= indicated by ROI block 105. After quantization, an inverse transform is applied to the remaining coefficients via 
55 transform block 106, followed by DC and optional component transform block 107. This results in generation of a 
reconstructed image. 

[0012] The JPEG2000 standard leaves many choices to implementers. 
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SUMMARY OF THE INVENTION 
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[001 31 A system comprising a memory and wavelet processing logic is descnbed. The memory .s sized o include 
ines to store a band of an image and additional lines. The wavelet processing logic compnses a wavelet trans orm 
Tnd access logic The wavelet transform generates coefficients when app.ied to data in the memory. The access logic 
iT^Liy into the line buffers to supply data stored in the memory to the wave.et transform and fo 
store coefficients in the memory, such that after data stored at a first pair of lines is read from memory into he buffers 
oftheaccesS 

that are associated with a second pair of lines different from the first pair of lines. 
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BRIEF nFSCRIPTION OF THE DRAWINGS 



[0014] The present invention will be understood more fully from the detailed description given I ; nd ^" the 
accompanying drawings of various embodiments of the invention, which, however, should not be taken to limit the 
75 invention to the specific embodiments, but are for explanation and understanding only. 
[0015] Figure 1 is a block diagram of the JPEG 2000 standard decoding scheme. 
r0016] Figure 2 illustrates one embodiment of an organization for an image in memory. 

0017 Figures 3A-F illustrate a transform memory organization for various levels depicting conceptually how coef- 
licients may be stored for the forward (Figures A-C) and inverse (Figures 3D-F) transforms. 
20 [0018] Figures 4A and B illustrate embodiments of a single memory where the input image data and the vanous 
decomposition levels of the image data can be stored during forward and inverse transforms, respectively. 
[0019] Figure 5 illustrates one embodiment of the process of handling the input data. 

[0020] Figure 6A illustrates a system having a progression order conversion parser. .^.....b,, 
[0021] Figu,e6Billustrate S aprogressionconverteroonvertingfromaresolutionprogress.vecodestreamtoaqual.ty 

25 progressive codestream. ^ . hor 

t0022] Figure 7A shows multiple ways to convert a codestream from one progression order to another 
[0023] Figure 7B shows one embodiment of simplified conversion paths to convert a codestream from one progres- 

sion order to another. _ w „ .„ 

r0024] Figure 8 illustrates one embodiment of a process for performing progression order convers.on. 
30 [o025] Figure 9 illustrates a decoder that selects portions of a codestream based on sideband information. 
[0026] Figure 10 is a flow diagram of a process for using layers when decoding. 
r00271 Fiqure 1 1 is a flow diagram of one embodiment of an editing process. 

[0028] Figure 12 illustrates a bell-shaped curve of a range of values that are quantized to a particular value. 
r00291 Figure 13 is a flow diagram of one embodiment of a process to reduce flicker. 
35 [0030] Figure 14 illustrates one embodiment of an encoder (or portion thereof) that performs the quantization to 
reduce flicker. 

[0031] Figure 15A illustrates a process for performing rate control. h„ m «« 
0032 Figure 15B illustrates an exemplary number of layers that may be subjected to first and second passes. 
r0033] Figure 16 illustrates one embodiment of the process for accessing the groupings of tile parts. 
[0034] Figure 17 and 18 illustrate quantizers for one component for a three level 5,3 transform. 
[0035] Figure 19 illustrates an example of HVS weighted quantization. 
[0036] Figure 20 is a block diagram of one embodiment of a computer system. 
[0037] Figure 21 illustrates an example progression with tile parts for a single server. 
[0038] Figure 22 illustrates an example of layering for a 5,3 irreversible transform. 

0039 Figure 23 illustrates an example in which transform has 5 levels and the data is divided up into layers 0-3. 
[So40] Figure 24 illustrates one example of a situation in which flicker may be avoided in which values in first and 
third trames are used to set the value in the second frame. 

[0041] Figure 25 is a block diagram of a prior art decoding process that includes color management. 
[0042] Figure 26 illustrates one embodiment of a non-preferred camera encoder. 
so [0043] Figure 27 illustrates one embodiment of a simpler camera encoder. 

[0044] Figure 28 is a flow diagram of one embodiment of a process for applying an inverse transform with clipping 
on partially transformed coefficients. 

DETAILED DESCRIPTION OF THE PRESENT INVENTION 

[0045] improvements to compression and decompression schemes are described. It is a purpose of «^"* u ~ 
and imputations described herein to use choices in JPEG 2000 lo make high speed, low cost, low memory and/ 
or feature rich implementations. 
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[0046] In the following description, numerous details are set forth in order to provide a thorough explanation of the 
present invention. It will be apparent, however, to one skilled in the art. that the present invent.on may be practiced 
without these specific details. In other instances, well-known structures and devices are shown in block diagram form, 
rather than in detail, in order to avoid obscuring the present invention. ^, mK ^i^ 
[0047] Some portions of the detailed descriptions which follow are presented in terms of algonthms and symbolic 
representations of operations on data bits within a computer memory. These algorithmic descriptions and representa- 
tions are the means used by those skilled in the data processing arts to most effectively convey the substance of their 
work to others skilled in the art. An algorithm is here, and generally, conceived to be a self -cons.sten sequence of 
steps leading to a desired result. The steps are those requiring physical manipulations of physical quantrt.es. Usually, 
though not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, trans- 
ferred combined, compared, and otherwise manipulated. It has proven convenient at times. pnncipaHy for reasons of 
common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like. 
[00481 It should be borne in mind, however, that all of these and similarterms are to be associated with the appropriate 
physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated othenv.se 
as apparent from the following discussion , it is appreciated that throughout the description, discuss.ons ut.l.z.ng temis 
such as "processing" or "computing" or "calculating" or "determining" or "displaying" or the like, refer to the action and 
processes of a computer system, or similar electronic computing device, that manipulates and transforms data repre- 
sented as physical (electronic) quantities within the computer system's registers and memories into other data s.milarly 
represented as physical quantities within the computer system memories or registers or other such information storage, 
20 transmission or display devices. . mav be 

[0049] The present invention also relates to apparatus for performing the operat.ons herein. This apparatus may be 
specially constructed for the required purposes, or it may comprise a general purpose computer selectively act.va ted 
or reconfigured bv a computer program stored in the computer. Such a computer program may be stored in a computer 
readable storage medium, such as, but is not limited to, any type of disk including floppy ^sks. optrca^d^s, CD-ROMs, 
and maanetic-ootical disks, read-only memories (ROMs), random access memones (RAMs). FPROM,. ^ pn ^ 
magnetic or optical cards,. or any type of media suitable for storing electronic instructions, and each coupled to a 

SrVhe 3 ?^ and displays presented herein are not inherenty related to any particular computer or other 
apparatus Various general purpose systems may be used with programs in accordance with the teachings here.n, or 
it may prove convenient to construct more specialiaed apparatus to perform the required method steps. The required 
structure for a variety of these systems will appear from the description below. In addition, the present invention is not 
described with reference to any particular programming language. It will be appreciated that a vanety of programming 
languages may be used to implement the teachings of the invention as descnbed herein .... „, h »m» 
[0051] A machine-readable medium includes any mechanism for storing ortransmitt.ng information in aform eadable 
by a machine (e.g., a computer). For example, a machine-readable medium includes read only memory ( ROM , 
random access memory ("RAM"); magnetic disk storage media; optical storage media; flash memory devices; electrical, 
optical acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital signals, etc.); etc. 
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Overview 

[0052] The following descriptions relate to implementations or novel ways to take advantage of the flexibility of JPEG 
2000 or other coding schemes with similar features. 

Memory Usage for Low Memory and Fast Burst Access 

r00531 Fiqure 2 shows one embodiment of an organization for an image in memory 201 Referring to Figure 2, only 
the "tile height" raster lines, or a band of the image, are in memory 201 , not the whole image. Thus, the amount of an 
image in memory 201 is equal to the image width multiplied by the tile height. Inside the band of the .mage .s at least 

WW Th C e h wavefeftransform processing logic 202 includes memory access logic 202A to read data from and store 
data to memory 201 to enable wavelet transform 202B to be applied to the data (image data or coefficients dependmg 
on the level of coefficient). Wavelet processing logic 202 may comprise hardware, software or a comb.nat.on of botlv 
f0055] In one embodiment, access logic 202A accesses the tile with four parameters: a pointer or index to the start 
of the tile in memory, the width of the tile, the height of the tile, and the line offset to get from the start of one line to 
another (the image width). Alternatively, access logic 202A accesses memory 201 using a pointer or index to the end 

of the tile instead of the width of the tile. 

[0056] In one embodiment, in order lo access for each line of a tile or a portion of a line of an image to perform some 
function F, the following process may be used. 
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line = start 

for y = 0 to tile_height - 1 

for x = 0 to tile_width - 1 

perform function F with line[x] 
line = line + line_offset 

One of the functions F may include applying a wavelet transform on pairs of lines. Also another function F may be a 

SST SS SSS by processing , g ic that may comprise hardware (e.g debtee Jljg, 

circuitry, etc.), software (such as is run on a general purpose computer system or a ded.cated machine), or a comb,- 

[o a 0 i 5 O 8 n ] O, «n 0 one embedment, coefficients from a subband are accessed using a similar scheme with a starting, **£ 
width height and line offset. Because rows of coefficients are stored together m memory, rows may be accessed 
eSnt I when stored in cache, burst accessible memories or memories that are wider than one coefficient. 
p059 Figures 3A-C show a transform memory organization for various .eve.s fP>^^^^^^ 
mav be stored All LH HL and HH coefficients (using the nominclalure of ITU-T Rec.T.800/ISO/IEC 154441 .2000 JPEG 
2000 ima e Cod ng System) are coded. These coefficients are represented by dotted lines in Figures 3B and C. Input 
Mnes of input tile 301 and LL coefficients (shown as solid lines in Figures 3B and 3C ,n successive e veto) only need 

be s ored temporarily while computing the transform w«h the exception of the final transform 
which are coded When a transform is used that does the horizontal and vertical transforms .n one pass and uses line 
buffers ^ once a pair of input rows has been complete* read (input lines or LL coefficients), the space used by lines 

" S Tg^es 3A-C show input tile 301 , level 1 (L1 )(302) and .eve. 2 (L2)(303) memory areas aligned with an off set 
o indicate how reuse might be accomplished in one embodiment. The addition of two rows, rows 312 and 313, to the 
m mon spec: used to hold input ti.e 301 , is needed to generate the L1 coefficients when -™^J^2lS. 
tile 301 for L1 coefficients. The addition of fourrows, rows 341-342, tothe memory space used to holdthe L1 
is needed to generate the L2 coefficients when reusing the memory storing the L1 coefficients for L2 coefficients^Note 
That there are two rows between rows 341 and 342 that are wasted space.) The additiona. .ines are preferab.y beh.nd 
the direction the wavelet transfoim is being applied to the information in the memory. 

[00611 More specifically, a pair of input rows input tile 301 corresponds to one row of each of LL LH HL and HH 

coeff ients at level 1 resulting from the application of a transform to two different l.nes and the results o applying the 

^STmoTm^ written into lines of the memory. For example, the results of applying a wavelet transform to 

Tnout rows sVo Jnd 311 are the coefficients in portions of rows 312 and 313 of L1 coefficients (302). ^xamp.e^LL 

cSSs 321 of row 312 corresponds to the LL coefficients (solid line) of level 1 , HL coeftoen* 322 of row 312 

cor^sponds fo the ^coefficients of level 1 . LH portion 323 of row 313 corresponds to the LH coefficents of level 1 , 

rndXort^n 324 Responds to the HH coefficients of level 1 . Note that the level 1 coefficients from the first two 

fnou Une ; a e stored in two extra rows at the top of the memory with the remaining level 1 coefficients being wntten 

£C£Z^Mg the data of input tile 301 to reuse the memory. The width and height for each type of coefficient 
nto the locauons to ng p ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

,0 the next LL row for level 1 (e.g., the distance from row 312 to row 314 in Figure 3B) is tw,ce the t.le w.dth (since 
parh coefficient row is from an area corresponding to two lines). 

4. SS^lSly. the results of applying the wavelet transform to two rows of LL coeffici 'eve. 1 

ire the coefficients in two rows namely LL coefficients (331). LH coefficients (332) HL«M ^eoffse, 
coefficients (334) at level 2. The width and height for level 2 coefficients is a quarter thai of input t.le 301 The line offset 
'or levSs?' r ttmes the ti.e width (since each coeffcient row is from an area corresponding to two level 1 LL rows 
or Z 1« lines). Thus, four extra lines of memory are needed to use the same memory that is .storing , the input t e 

so Z Ito e the 12 coefficients. Note that if a third decomposition level was being performed an add.t.ona 8 l.nes . would 
be needed Thus in this example, a total of 14 extra lines are needed to enable reuse of the memory that stores an 

TnpuT^ 

of extra lines is as follows: 
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[0063] To access subbands. such as the LL, LH, HL and HH subbands. only a starting pointer and the offset -between 
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rows/lines are necessary. The height and width are also needed to know when to stop when accessing a tile. 
[0064] As the number of decomposition levels increases, some rows at the bottom of memory become unused. That 
is. the lines of memory below the L1 coefficients after the first decomposition level become unused, the lines of memory 
below the L2 coefficients after the second decomposition level become unused, etc. In one embodiment, this extra 
space may be reused. 

[0065] Figures 3D-3F illustrate the corresponding inverse translorm memory usage in which additional lines store 
the results of applying an inverse transform and those additional lines are in the memory behind the direction the 
inverse transform is being performed. 

[0066] Figure 4A shows one embodiment of a single memory where the input and the various levels can be stored 
during application of a forward transform. Referring to Figure 4A, locations for the input tile, level 1 coefficients, level 
2 coefficients, and level 3 coefficients is shown with the added 2, 4 and 8 lines respectively. Figure 4B shows a similar 
single memory embodiment where the input coefficients of various levels of the transform can be stored along with the 
output during application of an inverse transform. 

[0067] Table 1 shows the amount of memory required for various transform levels for a 256x256 tile for separate 
memories and reused memory. 

Table 1 



level 


Separate memory (bytes) 


reused memory (bytes) 


1 


256x256 = 


65.536 


2x256 = 


512 


2 


128x128= 


16 ; 384 


4x256 = 


1,024 


3 


64x64= 


4,096 


8x256= 


2,048 


4 


32x32 = 


1 .024 


1 6x256 = 


4 ? 096 


5 


16x16= 


256 


32x256= 


8,192 




PvP - 

I 


64 


54x255 = 





[0068] For reused memory, the amount listed is the additional new memory used for that level. For this example, 
reusing memory for levels 1 , 2 and 3 saves memory. Level 4 may use a separate memory. 

[0069] The memory for levels 4, 5 and 6 could be placed in a single memory after level 3 has been generated or in 
a completely different and separate memory. The amount of memory necessary is 38x32, which is less than 5x256. 
Because there are two unused lines after generating the level 1 coefficients (i.e., the memory that stored the last two 
lines of input data), a small memory savings can be achieved by letting the levels 4, 5 and 6 reuse these two lines. 
This is particularly important because the number of additional lines for levels 4, 5, and 6 is 16, 32 and 64, and the 
extra space between the lines will be twice as far and half as wide as the level before. 

[0070] In one embodiment, coefficients from levels 4, 5, and 6 are packed in a smaller memory structure, such as 
storage area 450 in Figure 4. Referring to Figure 4, the level 4 coefficients are stored in an area having a height equal 
to the tile height divided by 8 (2^ where 3 corresponds to the number of levels) and a width equal to the tile width w 
divided by B (2 3 where 3 corresponds to the number of levels previously stored elsewhere). An additional two lines 
451 are all that is needed to store level 5 coefficients in the same necessary storage area. Similarly, an additional four 
lines is all that is necessary to accommodate using this memory storage area for the level 6 coefficients. Note that no 
lines are skipped when storing the coefficients. In one embodiment in which a 256x256 tile is being processed, the 
extra 5 lines at the bottom of storage area 430, two lines 421 and approximately 4.75 lines 422 are used to accommodate 
storage area 450. As shown, the approximate by three lines 422 represent allocated memory or in addition to that 
necessary to store the input tile. In this manner, the storage area for the input tile is almost completely reused. 
[0071] In one embodiment, to use a very little, or potentially minimum, memory, level 6 is stored separately from 
levels 4 and 5. However, this only saves 64 bytes of memory. 

[0072] A memory a little smaller than 273x256 can hold all the transform coefficients for a 256x256 tile. This is less 
than 7% more than a true in-place memory organization. Unlike an in-place memory organization, extra copies are 
avoided while simultaneously keeping the rows packed together for fast access. 

[0073] Table 2 shows another example of using separate versus reused memory for 128x128 tiles. For this size, the 
first three transform levels can reuse memory in a 142x128 buffer. 

Table 2 
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level 


Separate memory (bytes) 


reused memory (bytes) 


1 


128x128 = 


16 ; 384 


2x128= 


256 


2 


64x64= 


4.096 


4x128= 


512 
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level 


Separate memory (bytes) 


reused memory (bytes) 


3 


32x32= 


1.024 


8x128= 1024 



[0074] In one embodiment, a decision to use in-place memory or new memory is a function of tile height and transform 
level. Such a decision may be based on the following: 

if tile height > 2&* eveh2 K then use in-place method 
if tile height = 2^ tevel ' 2 > > then either may be used 
if tile height < 2( 3 * tevel_2 >, then use new memory 
15 To illustrate the application of the decision, Table 3 below: 

Table 3 



level 


2 A (3*level-2) 


1 


2 


2 


16 


3 


12B 


4 


1024 


5 


8192 



[0075] In some applications, adapting the memory organization to the tile height is inconvenient. A single fixed mem- 
ory organization can be used. Tile sizes smaller than 128x128 typically result in bad compression performance, so 
would typically not be used. While tile sizes bigger than 1Kx1K can be used for very large images, this does not 
significantly improve compression and the large amount of memory required would typically be burdensome. Therefore, 
assuming a tile height between 128 and 1024 inclusive and using in-place memory for 3 levels of the transform is a 
good heuristic. 

[0076] Decoding is similar in that the results of applying an inverse transform are written ahead of where the decoding 
processing logic is reading, with the only notable difference being that the start is from the highest level to the lowest 
level, such as level 6 to level 1 in the example above. In such a case, the input tile ends up at the top of the memory 
structure. The extra lines to accommodate the memory reuse are in decreasing order. For example, using the structure 
of Figure 4B, 8 lines would be necessary to create the L2 coefficients from the L3 coefficients, 4 extra lines would be 
necessary to create the L1 coefficients from the 12 coefficients and 2 extra lines would be necessary to create the 
input tile from the L1 coefficients. 

[0077] In one embodiment, to handle input tile data., a color conversion may be performed on the data prior to en- 
coding. Figure 5 illustrates one embodiment of the process of handling the input data. Referring to Figure 5, color input 
pixels are received in raster order. These color pixels may be in RGB. YCrCb, CMY, CMYK, grayscale, etc. The color 
input pixels may be stored as tiles in a memory, such as memory 501 , by band (or other forms). 
[0078] Pixels from storage 501 or received directly form the input undergo color conversion and/or level shifting, with 
the resulting outputs being stored in one coefficient buffers 502 r 502 N . That is, once the color conversion has been 
completed on each tile, it is stored in one of the coefficient buffers 502 r 502 N , and then the next tile can be processed. 
In one embodiment, there is one coefficient buffer for each component. 

[0079] Coefficient buffers 502 r 502 N are used by the transform in the manner described above to perform the wavelet 
transform while reusing memory. Thus, coefficient buffers 502 r 502 N are both input and output to wavelet transform. 
[0080] After the transform is applied to coefficient buffers 502 r 502 N , the context model 503 and entropy coder 505 
can perform further compression processing on the already transformed data. The coded data is buffered in coded 
data memory 505. 

[0081] While performing the further compression processing on one tile, the transform may be applied to another 
tile. Similarly, any or all the operations may be performed on multiple tiles at the same time. 
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Progression Order Conversion 

[0082] In the JPEG2000 standard, data in a compressed codestream can be stored in one of the five Progression 
orders. The progression order can change at different points in the codestream. The order is defined by embedded 
"for lavers" on layers, precincts, resolution, and components. 
[0083] Fiveprogressionorderca^ 

resolution-component-position progression (LRCP), resolution-layer-component-position progress.on (RLCP),i resolu- 
tion-position-component-layer progression (RPCL), position-component-resolution-layer progress.on (PCRL), compo- 
nent-position-resolution-layer progression (CPRL). W «f=..H 

[00841 The order may be defined in the COD or POC markers of the JPEG 2000 standard. The Coding style default 
/COD marker is defined by the JPEG 2000 standard and describes the coding style, number of decomposition evels, 
and layering that is the default used for compressing all components of an image (if in the main header) or a tile (if in 
a tile-part header). The Progression order change (POC) marker describes the bounds and progression ' °rderf°r any 
progression order other than that specified in the COD marker segments in the codestream. The Packet Length Main 
Header (PLM) indicates a list of packet lengths in tile-parts for every tile part in order and the Packet Length, Tile-part 
header (PLT) indicates tile packet lengths in a tile-part and indicates where the data is in the codestream^ 
r00851 The JPEG 2000 standard in section B.12 only specifies how packets of compress data are formedfor a given 
progression order It does not describe how data should be converted from one progression order to another progress.on 



[0086] In one embodiment, a progression order converting parser converts a codestream to a desired progression 
order based on the user input without decoding the data and then encoding it again. Figure 6A illustrates a system 
having such a parser. Referring to Figure 6A, parser 601 receives requests from a client for a particular progress.on 
order The client may be viewing a web page and selects a particular link. In response to the request, parser 601 
accesses server 602 to obtain the codestream associated with full image 603 from memory 604 and converts the 
codestream into a different progression order based on the request. The request indicates the progress.on order by 
using an optional command (e.g., RL2L (Resolution-layer progression to Layer Progression)). The progress.on order 
that is described may be based on layer, resolution, component, precinct, or tile. 

[0087] Figure 6B illustrates the progression converter converting from a layer progressive codestream (LRCP) to a 
resolution progressive (RLCP) codestream. The progression orders map directly to each other. 
r0088] Figure 7A shows multiple ways to convert a codestream from one progression order to another. Referring to 
Figure 7A. each of the five progressions (LRCP, RLCP, RPCL, CPRL, and PCRL) are shown with paths to each of he 
others such that all progressions are shown. In one embodiment, the parser causes all converses to go through the 
laver progression first and then to a selected conversion. Figure 7B shows one embodiment of such s.mpl.f.ed conver- 
sion paths in which the number of required mappings is reduced from 10 (as in Figure 7A) to 4. However, any one of 
the five progression orders could be used as the one to which all are converted before arriving at the selected order 
The conversion technique described herein simplifies source codes in that the number of lines of source code is much 
less than the multiple ways of conversion. This results in less debug time and fewer memory and run-t.me variables. 
[0089] To perform the conversion, the order of the packets in the codestream must be reordered^ The Packets are 
abeled by their sequential order in the codestream. Markers may indicate the starting point of the data, the length o 
the data (or alternatively the endpoint of the data) and how the data should be handled. For example, the .nd.cat.on of 
how the data is to be handled may indicate whether the data is to be deleted, whether the data is to be truncated, or 
some other operation to be performed on the data. Such handling information may also come from rale d.stortion 
intormation, such as may be provided in a PLT/PLM and/orthe PPT/PPM marker sets of the JPEG 2000standard. In 
this manner, the codestream may be truncated without changing the packet header. lnHir . atinn 
[0090] in one embodiment, a list, array, or other structure (such as reordering structure 601 A) is bu.lt by .nd.cat.ng 
the portion of data in each packet. Using this structure, the packets may be reordered. 

[00911 Figure 8 illustrates one embodiment of a process for performing progression order conversion. The process 
is performed by processing logic that may comprise hardware (e.g., dedicated logic, circuitry, etc.), software (such as 
is run by, for example, a general purpose computer or dedicated machine), or a combination of both. 
[0092] Referring to Figure B, the process begins by processing logic building a list from headers in the packets 
processing block 801 ) and optionally marking list items "delete" for quantization (processing block 802). Next, process- 
ing logic reorders the list to map the original progression to a desired progression (including handling input and output 
with progressions specified with POC markers (bounds on the progression order) (processing block 803). Thereafter, 
processing logic outputs coded data based on reordered list (processing block 804). 

[0093] Therefore, the combination of re-ordering and parsing allows specification of the desired order.ng and reso- 
lution, quality, etc. 
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A Progression Order Conversion Example 

[0094] The following is an example showing how packets are arranged in a codestream. The codestream was formed 
based on 2 components, 2 layers, 3 decomposition levels, and layer progression. 

[0095] Table 4 shows the packet order, length and association index of packets in the example. The packet order 
column shows the sequential order of packets placed in a codestream. The length indicates the length of the packets. 
The association index shows the resolution, layer, component and precinct of the packet. 

[0096] For example, packet[0] is the first packet in the codestream after the first tile header. It has a length of 589 
bytes. Association index RwLxCyPz indicates the packet belongs to resolution w, layer x, component y and precinct z. 

Table 4 



50 



Packet order 


Length 


Association Index 


packet[0] 


length=589 


ROLOCOPO 


packet[1] 


length=589 


R0L0C1 PO 


packet[2] 


length=924 


R1LOC0PO 


packet[3] 


length=924 


R1L0C1P0 


packet[4] 


length=1602 


R2LOC0PO 


packet[5] 


length=1602 


R2L0C1P0 


packet[6] 


length=733 


R3L0C0P0 


packet[7) 


length=733 


R3L0C0P0 


packet[8] 


length=535 


R0L1C0PO 


packet[9] 


length=535 


R0L1C1P0 


packet! 10] 


length=1523 


R1L1C0P0 


packet[11] 


length=1523 


R1L1C1P0 


packet[ 12] 


length=5422 


R2L1C0P0 


packet[ 13] 


length=5422 


R2L1C1P0 


packet! 14] 


length= 16468 


R3L1C0P0 


packet! 15] 


length=16468 


R3L1C1P0 



[0097] In this codestream, packets are grouped based on the layer in which they reside. The first 8 packets belong 
to Layer 0. The following 8 packets belong to Layer 1 . 

[0098] Using the conversion process described herein, the above codestream is converted to resolution layer pro- 
gression. The following shows how the above packets are re-ordered. 

[0099] After the layer progressive codestream is converted to resolution progression, in the new codestream, packets 
are grouped based on resolution. Such a grouping is shown in Table 5. The first 4 packets belong to resolution 0, the 
next 4 packets to resolution 1 , and so on. 



Table 5 
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Previous Packet order 


Packet order 


Length 


Association Index 


0 


packet[0] 


length=589 


ROLOCOPO 


1 


packet! 1) 


length=589 


R0LOC1PO 


8 


packet[2] 


length=535 


R0L1C0P0 


9 


packet[3] 


length=535 


R0L1C1P0 


2 


packet{4] 


length=924 


R1L0C0P0 


3 


packet(5] 


length=924 


R1L0C1P0 


10 


packet[6] 


lengths 1523 


R1L1C0PO 


11 


packet[7] 


lengths 1523 


R1L1C1P0 


4 


packetf8] 


length=1602 


R2L0C0P0 


5 


packet[9] 


length=1602 


R2L0C1P0 


12 


packet[10] 


length=5422 


R2L1C0PO 
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Tabic 5 (continued) 
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Previous Packet order 


Packet order 


Length 


Association Index 


13 


packet[11] 


iength=5422 


R2L1C1P0 


6 


packet[12] 


length=733 


R3L0C0P0 


7 


packet[13] 


length=733 


R3L0C1 PO 


14 


packetf14] 


length=16468 


R3L1C0P0 


15 


packet[15] 


length=16468 


R3L1C1P0 



20 



One Embodiment of a Conversion Algorithm 
[01 00] Resolution to Layer Progression 

n = 0; 

for(l=0;l<layer;l++){ 

for(r=0;r<resolution+l;r++){ 
for(c==0;c<component;c++){ 

new_packet[n] = old_packet[l*component + r*iayer*component + 
c]; n++; 
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Layer to Resolution Progression 

n = 0; 

for(r=0;r<resolution+l;r++){ 
for(l=0;l<layer;l++){ 

for(c=0;c<component;c++) { 
new_packet[n] = old_packet[r*component + 

l*(resolution+l)*component + c]; 
n++; 

) 

} 



55 



} 

where layer = the number of layers in a codestream, 

resolution = the number of decomposition levels in a codestream, and 

component = the number of components in a codestream 

Data Hiding (Sideband Information) in JPEG2000 Coding 

[0101] Bit hiding allows sideband information to be transmitted without increasing the file size. Sideband information 
that does increase file size but does not break naive decoders might also be valuable (although the COM marker 
defined by the JPEG 2000 standard might be used instead). 
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[01021 Some marker segments, packet headers and packets are padded out to the nearest byte. Examples of the 
JPEG 2000 marker segments include PPM, PPT, PLM. and PLT. In addition, some marker segments can be longer 
than they need to be including QCD, QCC, and POC. In all of these cases, the padded data values are not defined. 
[01031 Several proprietary coding schemescould use this semi-randomly located undefined data to prov.de a number 
5 of important types of information including, but not limited to, decoding and filtering hints, ownership segmentation 
hints and so on. A hint might include an index to a particular enhancement scheme. For example, rf rt is > known hat 
an image is mostly text, a value may be sent that indicates that a first post-processing filter is to be used. On the other 
hand, if the area is mostly a graphic image, then a value may be sent that indicates that a second post-processing 

io [0104] t0 Th e e"ollowing are places where bits may be hidden or sideband information may be stored in the codestream. 

arithmetic coder (AC) termination (without predictable termination) 
end of packet header rounding to byte 
after last packet, before next tile 
tag tree construction by not always using minimum 
packet header Lblock signalling 

LSB parity for codeblocks (refinement pass only, cleanup pass only, all) 
QCD, QCC extra subbands, POC. 

[0105] For example, with respect to hiding data using AC termination, 0 to 7 bits are provided, at least everytirne 
the coder is terminated. However, this could be extended for a few bytes. These extra bits and bytes may be used for 

sending extra information. Tkft - Arfl 

r0106] With respect to each packet header, the end of a packet header is rounded to a byte boundary. Therefore, 
there may be 1 to 7 bits that may be available for sending extra information at times when rounding would have been 
necessary. Similarly, each packet is rounded to a byte boundary, thereby providing 1 to 7 bits (assuming that rounding 
would have been necessary). Also the last packet in a tile-part can be extended a few bytes. These extra bytes may 
be used to send additional information. 

[01 07] The length of the compressed data for a code-block can be given in the packet header with a non-minimum 
representation. The choice of representation (e.g., a non-minimum representation) could be used for indicating other 

[0?Ml 8tl with respect to tag tree data hiding, packet headers of the JPEG 2000 standard use tag trees for coding first 
inclusion and zero bitplane information. When there are multiple codeblocks, tag trees are like a quadtree of minimum 
values For example, in the case of 1 6 codeblocks in a 4x4 arrangement in a packet, the arrangement may be as follows: 
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10 


7 


12 


15 


3 


20 


21 


5 


81 


45 


5 


9 


18 


8 


12 


24 



An example tag tree, which is minimal for the 4x4 arrangement above is as follows: 



3 


0 


2 


7 


4 


7 


10 




5 


2 


0 


17 


16 


0 








73 


37 


0 


4 








10 


0 


7 


19 



in which "3" is added to every codeblocks value, and »0», "2», "5» and "2" are each added to the 4 « respon ding 
codeblocks. Finally, there is one value per codeblock. That is, the minimal tag tree is created by taking the firet 2x2 
group in the 4x4 arrangement above and look at minimum value is out of the four values. In this case, for the 2x2 block 



10 


7 


3 


20 



the minimum value is 3. This is then performed on the other 2x2 blocks. Then these identified minimum values are 
evaluated again to determine their minimum, which would be "3" in the example. Then the minimum value is subtracted 
from the four minimum values to create the following 
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0 
5 



2 
2 



Then for the remaining numbers in the 4x4, the number 3 is subtracted from each value along with the value in the 
2x2 that corresponds to the particular value in the 4x4 arrangement, thereby resulting in the tag tree above. 
The first row adds up as follows: 



10 



10=3+0+7 



7=3+0+4 



15 



12=3+2+7 



15=3+2+10 
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[0109] A variable length code may be used that efficiently represents small numbers. 
[0110] An example of a tag tree that is not minimal is as follows: 



1 3 

63 



7 
u 

73 
10 



4 
1/ 
37 

0 



7 

16 
0 
7 



10 
0 
4 

19 



(Note that representing "3", "0", "2", "5" and -2" might use less bitstream data than "2", "1 V3", "6" and "3 



r01 1 1 ] Once a tag tree representation has been made, a determination can be made as to whether the representation 
is minimal or not based on whether there is a zero in the 2x2 block. Therefore, this information is hidden. For example, 
the 1 bit block represents the 1 in the 2x2 block above indicates it is not part of a minimal tag tree, but can be used to 
convey some particular information to a decoder. Likewise if a 2 was the minimal value in the 2x2 block, such a fact 
may convey different information to a decoder. 

[01121 The JPEG 2000 POC, QCD, and QCC markers can have redundant entries. It is as if the codestream were 
quantized and the markers were not rewritten. For example, the QCD and QCC markers have values for a number of 
subbands specified by the syntax of the marker. If there are fewer subbands actually coded in the brtstream data may 
be hidden in the values used for the missing subbands. The redundant entries may be replaced and used for hidden 

or sideband information. . . 

[0113] The hidden or sideband information may include post-processing hints (such as, for example, sharpen this 
tile with a specified filter or strength, or smooth, or perform optical character recognition (OCR) on this region, etc.), 
decoding hints, security (such as, for example, an encryption key for decoding the remainder of the .mage or another 
image, etc.) codestream identification (such as, for example, labeling POTUS as the originator of the file, etc.) and/or 
other information. 

Use of Layers When Encoding 

[0114] Layers are part of the JPEG standard. In one embodiment, sideband information, possibly in a COM marker, 
is used by the decoder to allow selecting of layers during decoding. The sideband information may be used to select 
layers for postcompression quantization to meet rate/distortion targets for different viewing distances, different reso- 
lutions different regions of interest, different frequency content for analysis (e.g., finding edges of text). 
[0115]' in one embodiment, the layers are predefined based on rate. For example, the first layer represents a 1-b.t 
per pixel image, while the second layer represents a 2-bit per pixel image, etc. Therefore, the layers run from the lowest 
quality to the highest quality. Likewise, target rates can be met for lower resolutions as well. 

[0116] The sideband information may be stored in a marker segment of the codestream. In one embodiment, the 
JPEG 2000 comment (COM) marker is used to provide information about the layers. Specifically, the COM marker may 
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be used to indicate the number of bytes tor each resolution and/or rate across the entire image or a relative number 
of b"es layer. Table 6 indicates each layer and rts resolution in the number of bytes across the ti.e 

in an image. Such stable may have distortion values instead. 

Table 6 
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lev=0 


layer=0 


comp=0 


bytes=529 


lev=0 


layer=0 


comp=1 


bytes=555 


lev=0 


layer=0 


comp=2 


bytes=493 


lev=0 


layer=1 


comp=0 


bytes=129 


lev=0 


layer=1 


comp=1 


bytes=130 


lev=0 


layer=1 


comp=2 


bytes=123 


lev=0 


layer=2 


comp=0 


bytes=7 


lev=0 


layer=2 


comp=1 


bytes=8 


lev=0 


layer=2 


comp=2 


bytes=12 


lev=0 


layer=3 


comp=0 


bytes=1 


!ev=0 


layer=3 


comp=1 


bytes=1 


!ev=0 


layer=3 I 


comp=2 


bytes=129 


I lev=1 


layer=0 


comp=0 


bytes=705 


lev=1 


layer=0 


comp=1 


bytes=898 


lev=1 


layer=0 


comp=2 


bytes=712 


lev=1 


layer=1 


comp=0 


bytes=146 


lev=1 


layer=1 


comp=1 


bytes=114 


lev=1 


layer=1 ] 


comp=2 


bytes=116 


lev=l 


layer=2 


comp=0 


bytes=224 


lev=1 


layer=2 


comp=1 


bytes=250 


lev=1 


layer=2 


comp=2 


bytes=263 


lev=1 


layer=3 


comp=0 


bytes=201 


lev=1 


layer=3 


comp=1 


bytes=212 


I lev=1 


layer=3 


comp=2 


bytes=200 


I lev=2 


layer=0 


comp=0 


bytes=889 


lev=2 


layer=0 


compel 


bytes=1332 


lev=2 


layer=0 


comp=2 


bytes=1048 


lev=2 


layer=1 


comp=0 


bytes=240 


lev=2 


layer=1 


comp=1 


bytes=329 


lev=2 


layer=1 


comp=2 


bytes=328 


lev=2 


layer=2 


j comp=0 


bytes=599 


lev=2 


layer=2 


comp=1 


bytes=767 


lev=2 


layer=2 


comp=2 


bytes=725 


lev=2 


layer=3 


I * _ 
comp=0 


bytes=335 


lev=2 


layer=3 


comp=1 


bytes=396 


lev=2 


layer=3 


I mmD=2 


bytes=420 


lev=3 


layer=0 


! comp=0 


bytes=1 


lev=3 


layer=0 


comp=1 


bytes=395 


lev=3 


layer=0 


comp=2 


bytes=402 


lev=3 


layer=1 


comp=0 


bytes=251 


lev=3 


layer=l 


comp=1 


bytes=450 


lev=3 


layer=1 


comp=2 


bytes=562 


lev=3 


layer=2 


comp=0 


bytes=525 


lev=3 


layer=2 


comp=1 


bytes=990 


lev=3 


layer=2 


comp=2 


bytes=1313 


| lev=3 


layer=3 


j comp=0 


bytes=1214 
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Table 6 (continued) 



lev=3 


layer-3 


comp-1 


bytes- 1798 


iev=3 


layer=3 


comp=2 


bytes=2585 
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[0117] In another embodiment, the ordering could be by layer. Thus, the information above is consolidated foreach 
level (not segregated by level or component), as shown below: 

Ordering by layer=0 bytes=7959 bitrate=0.971558 PSNR=30.7785 
Ordering by layer=1 bytes=10877 bitrate=1 .327759 PSNR=32.0779 
Ordering by layer=2 by tes= 16560 bitrate=2.021484 PSNR=35.7321 

[0118] Distortion by layers can be based on PSNR. For example, 

layer-0 PSNR-30.7785 
layer=1 PSNR-32.0779 
laycr=2 PSNR=35.7321 

[0119] In an alternative embodiment, such information may be hidden in the codestream as described above. The 
information may be used to control rate distortion. 

[0120] In another embodiment, the layers may be predefined for a particular viewing distance. In such a case, the 
data is divided into layers from the highest frequency, lowest resolution to the lowest frequency, highest resolution. 
[0121] In one embodiment, the layer information indicates the summation of bits across the entire image for that 
layer and all previous layers (for example the 16 f 011 bits listed next to layer 1 indicates the total number of bits for 
layer 0 and layer 1). Alternatively, bytes, words, kilobytes, or other units of memory or rate could be used instead of 
bits. Table 7 shows this type of absolute rate information. 

[0122] Table 8 shows relative rate information. Layer 0 has 4096 bits, layer 1 has 11 ,915 bits, etc. 

Table 7 
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layer 


Rate (bytes) 


0 


4,096 


1 


16,011 


2 


40,000 


3 


100,000 


4 


250,000 


5 


500,000 


6 


1 ,000,000 


7 


2,500,000 


8 


5,500,000 


Table 8 


layer 


Rate (bytes) 


0 


4,096 


1 


11,915 


■ 2 


23,989 


3 


60,000 


4 


150,000 


5 


250,000 


6 


500,000 


7 


1 ,500,000 


8 


3 r 000.000 



[0123] For example, if only 750,000 bytes may be allowed in the decoded image, then all that can be decoded (as 
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the 1,000,000 bytes tabulated with layer 6 includes the 500,000 bytes of layers 0-5) is through layer 5 and half of 
importance layer 6. In some embodiments, no packets from layer 6 would be included. In other embodiments, some 
packets from layer 6 would be included and others would be replaced by zero packets so that the total amount of layer 
6 data was approximately 250,000 bytes. 
5 [0124] Figure 22 illustrates an example of layering for a 5,3 irreversible transform with three levels, MSE or similar. 
Referring to Figure 22, there are 45 layers shown. Each additional layer improves MSE in an order that gives good 
rate-distortion for MSE. 

[0125] Figure 23 illustrates another example in which transform has 5 levels and the data is divided up into layers 
0-3. Layer 0 corresponds to the thumbnail version, layers 0-1 correspond to the monitor (or screen) resolution, layers 
10 0-2 correspond to the print resolution, and layers 0-3 correspond to lossless. 

[0126] In an alternative embodiment, the layers may be predefined for some other distortion metric (e.g., MSE, 
weighted MSE, sharpness of text, etc.) 

[0127] The decoder uses the information regarding the layers from the codestream to select layers to generate an 
image. The decoder knowing what the desired viewing characteristics from the application or implementation (see 

'5 Table 9 below), and using the information from the codestream specifying the layers, can quantize the codestream in 
order to display an image at the correct viewing distance. Figure 9 illustrates such a decoder. Referring to Figure 9, 
decoder 901 receives a codestream and includes quantization logic 902 that examines the COM marker and uses 
information about the viewing distance it is at stored in storage 903 to generate quantized codestream 904 via, for 
example, selecting the proper layers. Quantized codestream 904 is decoded by decoding logic 905 (e.g., a JPEG 2000 

20 decoder) after selecting layers to generate an image data 906. A naive decoder would simply ignore the data in the 
comment marker. 

[01 28] Figure 1 0 is a flow diagram of a process for using layers when decoding. The process is performed by process- 
ing logic that may comprise hardware (e.g., dedicated logic, circuitry, etc.), software (such as is run by, for example, a 
general purpose computer or a dedicated machine), or a combination of both. 
25 [0129] Referring to Figure 10, the process begins by processing logic receiving a codestream of -compressed logic 
data (processing block 1001). The image data is organized into multiple layers, each of which comprises coded data 
that adds visual value to the image (e.g., look sharper, better defined, better contrast, etc.). Next processing logic 
selects one or more layers for quantization based on sideband information (processing block 1002). After selection, 
processing logic decompresses the non-quantized layers of the codestream (processing block 1003). 

30 

Editing of Tiles, Tile-parts, or Packets 



[0130] Once a' codestream is created, it may be desirable to edit parts of the image. That is, for example, after 
performing encoding to create the codestream, a set of tiles may be decoded. After decoding the set of times, editing 

35 may be performed, followed by encoding the set of tiles with the edits to the same size as the encoded tiles were prior 
to their decoding. Examples of typical editing include sharpening of text and removing "red-eye." The JPEG 2000 
codestream can be edited in memory or in a disk file system without rewriting the entire codestream. 
[0131] Figure 11 is a flow diagram of one embodiment of an editing process. The process is performed by process 
logic that may comprise hardware (e.g., dedicated logic, circuitry, etc.), software (such as is run by, for example, a 

40 general purpose computer or a dedicated machine), or a combination of both. 

[0132] Referring to Figure 11 , processing logic initially determines the tiles, tile-parts, or packets that coverthe area, 
resolution, components, and/or precincts to be edited and decodes them (processing block 1101). This determination 
may be made in response to a user selecting an area and/or working resolution. The determination may use editing 
information for a higher resolution to determine which parts or tiles coverthe portion to be edited. Once decoding has 

45 been completed, processing logic performs the desired edits (processing block 1102).. 

[0133] After performing the desired edits, processing logic recompresses the data into coded data (processing block 
1103) and creates a replacement tile, lile-parl, or packet for the codestream (processing block 1104). In one embodi- 
ment, in creating the replacement tile, tile-part, or packet, processing logic pads out the data with bytes at the end of 
the codestream if the new data is smaller than the unedited version of the data to make the replacement tile, tile-part 

50 or packet the same size as the unedited version. 

[0134] In an alternative embodiment, processing logic may use a marker, or tag, such as a COM marker segment 
of the appropriate length instead of the padding. The COM marker could be used to fill space or could contain information 
that the encoder wanted to include. It could contain information such as, for example, sideband information described 
herein or a copyright license for an image or text or other file format information. 

55 [0135] In one embodiment, in creating the replacement tile, tile-part, or packet, processing logic truncates the last 
packets for any or all components until the data fits in the codestream if the new data is larger than the unedited version 
of the data. 

[0136] Editing of an image may be performed by changing coded data for tiles, tile-parts, or codeb locks. In one 
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embodiment, editing is performed without changing file size by quantizing instead of expanding. In another embodiment, 
a predetermined amount of extra space is allocated per tile or per codeblock to allow for a predetermined amount of 
expansion. In still another embodiment, coded data may be put at end of files by manipulating tile headers and putting 
invalid tile data in COM markers. 

5 [0137] Note that if there are subsequent tile-parts that depend on the data in the portion of the codestream that is 
being edited, these tile-parts may become useless in the codestream. An indication of this useless data may be noted 
to the decoder by one of several methods. These methods involve inserting or modifying information in the codestream 
to indicate the presence and/or location of the useless data. In one embodiment, the application uses a status buffer 
to indicate that the data in tile-parts subsequent to an edited tile-part may be useless. The status buffer may be in 

io workspace memory and describes dependencies between packets. If an earlier packet is altered, the subsequent 
packets cannot be decoded as is. These subsequent packets must be edited accordingly or eliminated. In another 
embodiment, such an indication may be made by zeroing out the data section of those tile-parts and/or creating a PPT 
marker segment that denotes no data. 

*5 Optimal Encoder Quantization 

[0138] During encoding, unquantized coefficients from some or all subbands maybe divided by a value of Q to create 
the quantized coefficient values. This value Q may have a wide range of values. Typical encoders quantize a number 
of the values in a single particular range of values is made equal to one single coefficient value. In essence, all the 

20 coefficients in the particular range are quantized to the same value. This can be exemplified by Figure 1 2 which shows 
that the range of values is often in a bell shaped curve and that all of the values in the particular range, such as range 
fl 1 are sent to the decoder as one quantized value, such as R , and the decoder will reconstruct these values to a 
particular value. Assume a decoder reconstructs these values to a predetermined value (e.g., floor ("& min + 1 /£ max), 
or min + Q, where Q is the quantization step size). For example, if the range of values is between 16 and 31 , then 

25 the decoder may assume the value is 24. In one embodiment, instead of using 1/2 as the value, another value is 
selected, such as floor (3/8 min + 5/8 max), or min + 3/8Q, where Q is the quantization step size. Therefore, if the 
range is from 1 6 to 31 , then it is assumed that the decoder will reconstruct the value to 22, instead of 24. 
[0139] In some cases, two spatially adjacent coefficients may be close to each other numerically yet in separate 
quantization bins, such as coefficient values 1201 of range R 2 and 1202 of range R n in Figure 12. The results of the 

30 quantization may cause an artifact to occur. In one embodiment, for coefficients near a boundary between two quan- 
tization bins, the encoder selects a bin such as Range R., into which a coefficient, such as coefficient 1201, will be 
quantized so that it is consistent with neighbors, such as coefficient 1202. This helps avoid artifacts. That is, this 
technique reduces distortion yet may increase rate, particularly when a coefficient is moved from a smaller bin to a 
higher bin. 

35 

Flicker Reduction for Motion JPEG 

[0140] At times, flicker occurs when applying wavelet compression to motion sequences. An example of such flicker 
may include the image getting brighter or darker in areas or the appearance of edges changing in successive frames 
40 as the motion sequence is played (mosquito noise around the edges). The flicker may be due to the application of 
different local quantization to successive frames of a motion sequence or to noise exacerbated by quantization that is 
viewed temporarily. 

[0141] To reduce flicker, coefficients that are in the same position and close to the same value in successive frames 
are forced to the same value. That is, the coefficients values in successive frames are set to a predetermined value. 
This is essentially a form of quantization that is applied during encoding. Figure 13 is a flow diagram of one embodiment 
of a process to reduce flicker. 

[0142] A test of whether to apply such quantization to a coefficient value in a subsequent frame is based on the 
quantization that was performed on the coefficient in the previous frame. Thus, the encoder is utilizing frame depend- 
ency to eliminate flicker while the decoder decodes data frame by frame independently. 
50 [0143] In one embodiment, in order to reduce flicker in motion JPEG, coefficient values are modified (quantized) 
based on their relationship with each other with respect to a threshold. For example, if Dn and Dn+1 are the corre- 
sponding coefficient (same spatial location and same subband) in two frames before quantization, if D'n and D'n+1 
represent these coefficients after quantization, if Q(») are scalar quantization, and if the value T is a threshold, then 
the following may be applied: 

55 

if (IQ(Dn+1)-(D'n) I < T) 



16 



EP 1 233 624 A1 



D'n+1 = D'n 



else 

5 

D'n+1 = Q(Dn+1) 

For example, the value T may be twice the quantization step size. Other values of T include, but are not limited to, 
10 J2Q : 1.50.272 0. 

[0144] One of the coefficient values may be modified to be either a predetermined closeness to another coefficient 
value. The closeness may be determined by some threshold. The threshold may be user set or adaptive based on 
some criteria. The threshold could be different based on the subband and, perhaps, on the persistance of the particular 
value (number of frames that this coefficient is close). In one embodiment, the coefficient value is set equal to the other 
15 coefficient value. In alternative embodiments, the coefficient is set to be within the quantization bin size of the other 
coefficient value or twice the quantization bin size. 

[0145] Figure 14 illustrates one embodiment of an encoder (or portion thereof) that performs the quantization de- 
scribed above. Referring to Figure 14, a quantizer 1400 receives coefficients 1410 for frames of a motion sequence 
from a wavelet transform (not shown). The coefficients are received by quantization logic 1401 which compares a 
20 threshold value stored in memory 1401 to coefficient values for the previous frame that are stored in memory 1403 to 
coefficients 1410 with a scalar quantizer Q applied from memory 1404. 

[0146] Quantization logic 1401 may comprise comparison hardware (e.g., logic with gates, circuitry, etc.) or software 
to perform the comparison. This comparison hardware and software may implement a subtractor or subtraction oper- 
ation. The results are a quantized codesteam (assuming some values have been changed.) 
25 [0147] This may be applied over two or more frames. Also the comparison is not limited to two consecutive frames. 
The comparison can be over 3, 4, 5, etc., frames, for example, to determine if a variance exists. Figure 24 illustrates 
one example in which values in a first and third frame are used to set the value in the second frame. 
[0148] Note that the quantization can also be codestream quantization with a code block-based rule. 

30 Rate Control, Quantization, and Layering 

[0149] In one embodiment, selective quantization of coefficients can be performed during encoding by setting a 
subset of the refinement bits to be the more probable symbol (MPS). This may be performed at a user selected bitplane. 
For examples, if there is text on a background image, with a goal of having sharp text images while minimizing coded 
35 data required for the background, the refinement bits that are set to MPS are those that do not effect text for the last 
bitplane, while using the actual value for bits that effect text. 

[0150] Such a quantization scheme may be used to implement non-uniform quantization step sizes. For example, if 
one wanted to have a background with fewer bits, setting the refinement bits to the MPS could operate as a form of 
quantization. This quantization scheme causes some level of distortion but lowers the bit rate necessary to transfer 
40 the codestream. 

[0151 ] Note that although this technique may be applied to bits generated during the refinement pass, the technique 
has application to other compression schemes (e.g., lists generated during subordinate passes, tail bits of CREW of 
Ricoh Silicon Valley, Menlo Park, California, MPEG IV texture mode, etc.). 

[0152] In one embodiment, the same technique may be applied to other changes between frames. That is, in one 
45 embodiment, a change due to a rate distortion in one frame may be performed in a subsequent frame to avoid distortion 
effects. 

Rate Control and Quantization 

50 [0153] In one embodiment, user specified quantization is provided. For a 3 level transform for one component, 7 
quantization values are sufficient: level 1 HH, level 1 HL and LH, level 2 HH, level 2 HL and LH, level 3 HH, level 3 HL 
and LH, and level 3 LH. 

[0154] If quantization values are bitplanes to truncate (which is equivalent to scalar quantization by powers of 2), 
3-bit values (0...7) are sufficient for most applications. (For image components with depth 12-bits or more and 5 or 
55 more transform levels, perhaps higher quantizations might be useful.) Values 0...6 could be used to specify the number 
of bitplanes to truncate and 7 could be used to mean discard all bitplanes. The three bit values may be written to a 
controller that controls compression (or decompression) hardware (e.g., JPEG2000 compatible hardware) to perform 
the quantization. 
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[0155] For 3 component color quantization: 

• 21 values can be used with separate values for each component, 

• 14 values can be used, 7 for luminance and 7 for chrominance, 

5 • 17 values can be used for 4:1 :1 subsampled data, 7 for luminance and 5 for each chrominance component, 

• 12 values can be used for 4:1 :1 subsampled data, 7 for luminance and 5 for chrominance, 

• 19 values can be used for 4:2:2 subsampled data, 7 for luminance and 6 for each chrominance component, and 

• 13 values can be used for 4:2:2 subsampled data, 7 for luminance and 6 for chrominance. 

10 Since 21*3= 63 bits is less than 8 bytes, transferring or storing the quantization uses little resources. A central process- 
ing unit (CPU) might select one predetermined quantizer from a table and write it to a CPU or other controller controlling 
special purpose JPEG 2000 hardware (a chip) for each frame of a motion JPEG 2000 video sequence. Alternatively, 
one implementation of JPEG 2000 might have a small memory that holds 8 or 16 different quantizers that could be 
selected for each frame. 

J5 [0156] Quantizers can also be used to assign bitplanes to layers. For example, Q 0 , Q 1( and Q 2 may be quantizers 
that specify bitplanes of coding pass to quantize. Quantizer Q 0 causes the most loss, while quantizer Q 2 causes the 
least loss. Layer 1 is ail the data quantized by Q 0 but not quantized by Q v Layer 2 is all the data quantized by Q., but 
not quantized by Q 2 . Layer 3 is all the data quantized by Q 2 . 

20 Simple Quantization 

[0157] Figures 17 and 18 show example quantizers (label A...Q) for the 3-level 5/3 transform as the number of 
coefficient LSBs to truncate or not code. Truncating N bitplanes is equivalent to a scalar quantizer of 2 N . The subband 
where the quantization changes with respect to the previous quantizer is highlighted with a dashed box. The quantizers 
25 D. K and Q all have the same relationship between the subbands. Other quantizers might be used that are better for 
MSE or for other distortion metrics, 

[0158] The exemplary Verilog below converts a single quantization value "q" into seven quantizers (number of LSBs 
to truncate). The variable q_1_HH is used for level 1 HH coefficients, the variable q_1_H is used for level 1 HL and 
LH coefficients, etc. Some consecutive values of q result in the same quantizer: 0 and 1; 2 and 3; 4 and 5; 8U-6 and 
30 8i+7 for all integers i with i > 0. 
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module makeQ(q, q_lHH, q_lH, q_2HH,-q_2H, q_3HH,q 3H, 
q_3LL); 

input [5:0] q; 
output [3:0] q_lHH; 
output [3:0] q_lH; 
output [3:0] q_2HH; 
output [2:0] q_2H; 
output [2:0] q_3HH; 
output [2:0] q_3H; 
output [2:0] q_3LL; 

wire [3:0] temp_2H; 
wire [3:0] temp_3HH; 
wire [3:0] temp_3H; 
wire [3:0] temp_3LL; 
wire [2:0] qlo; 
wire [2:0] qhi; 

assign qlo = q[2:0]; 
assign qhi = q[5:3]; 

assign q_lHH = qhi + ((qlo >= 2) ? 1 : 0); 
assign qJLH = qhi + ((qlo >= 4) ? 1 : 0); 

assign q_2HH = qhi + ((qlo >= 6) ? 1 : 0); 
assign temp_2H = qhi + ((qlo >= 1) ? 0: -1); 
assign temp_3HH = qhi + ((qlo >= 3) ? 0: -1); 
assign temp_3H = qhi + ((qlo >= 5) ? 0: -1); 
assign temp_3LL = qhi - 1 

assign q_2H = (temp_2H < 0) ? 0 : temp_2H; 
assign q_j3HH =' (temp_3HH < 0) ? 0 : temp_3HH; 
assign q_3H = (temp_3H < 0) ? 0 : temp_3H; 
assign qJ3LL = (temp_3LL < 0) ? 0 : temp_3LL; 



endmodule 



Human Visual System Weighting for Color and Frequency 



[01 59] Table 9 shows additional bitplanes to quantize (e.g. , truncate) for luminance to take advantage of the frequency 
response of the Human Visual System (from Table J-2 of the JPEG 2000 standard). A viewing distance of 1000 pixels 
50 might be appropriate for viewing images on a computer monitor. Larger viewing distances might be appropriate for 
print images or television. 



55 



3NSD0CID: <EP 1233624 A 1J_> 



19 



EP 1 233 624 A1 



Table 9 - 



5 



Human Visual System Weighting for Luminance 


subband 


extra biplanes to quantize for viewing distance of ... 


1000 pixels 


2000 pixels 


4000 pixels 


1HH 


2 


4 or 5 


discard all 


1HL, 1LH 


1 


2 or 3 


6 


2HH 




2 


4 or 5 


2HL, 2LH 




1 


2 or 3 


3HH 






2 


3HL, 3LH 






1 



75 Additionally chrominance may be quantized more heavily than luminance. 

[0160] Figure 19 shows a quantization that starts with Figure 17(D) and then adds frequency weighting for a 1000 
pixel viewing distance (to both luminance and chrominance), keeps 3LL chrominance unchanged, discards 1HL and 
1HH chrominance for 4:2:2 and additional 2 bitplanes are discarded for the remaining chrominance. 
[0161] Sharp text without ringing artifacts is more desirable than exact gray value for text/background. That is, if a 

20 gray level is supposed to be at 50% (for example), and is instead at 60%, it is often not visually objectionable if the 
image is of text. In one embodiment, the LL (DC) coefficients are quantized more heavily for text than for non-text 
images at low bitrate. For example, tor an 8-bit image component, a quantiation step size of 8, 1 6 or 32 might be used 
for text only regions and a quantization step size of 1 , 2 or 4 might be used for regions containing non-text. This allows 
more fidelity for the high frequency coefficients, thereby resulting in text with sharp edges. 

25 

Using Quantizers to Divide Things into Layers 

[0162] Table 1 0 shows 1 6 example quantizers. Quantizer 1 5 is lossless. Quantizer 8 is the same as Figure 19. These 
can be used divide the subband bitplanes into layers. 

30 

Table 10 





subband 


0 


1 


2 


3 


4^ 


5 


6 


7. 


8 


9 


10 


11 


12 


13 


14 


15 




Y 1HH 


ail 


all 


6 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


35 


Y 1HL, 
LH 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 




Y2HH 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 


0 


0 




Y2HL, 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 




LH 


































40 


Y3HH 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 




Y3HL ; 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 




LH 




































Y3LL 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


«0 


0 


0 


45 


C, 1HL, 










HL and HH always discarded for 4:1 :1 or 4:2:2 only 










HH 




































C 1 1 LH 


all 


ail 


all 


all 


6 


6 


; 5 


5 


4 


4 


3 


3 


2 


2 


1 


0 




2HH 


ail 


6 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 




2HL, 


all 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


50 


LH 




































3HH 


all 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 




3HL, 


all 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 




LH 
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C 1 3LL 


all 


0 


0 


o 


0 


0 


0 


0 


0 


0 


O 


0 


0 


0 


0 


0 


C 2 1HL. 










HL and HH always discarded for 4:1 :1 or 4:2:2 only 










HH 
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Tabic 10 (continued) 



subband 


0 


1 


2 


3 


4 


5 


6 


7. 


8 


9 


10 


11 


12 


13 


14 


1 ^ 


C 2 1LH 


all 


all 


ail 


all 


6 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


0 


C 2 2HH 


all 


6 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


C 2 2HL, 


all 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


LH 


































C 2 3HH 


all 


6 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


C 2 3HL, 


all 


5 


5 


4 


4 


3 


3 


2 


2 


1 


1 


0 


0 


0 


0 


0 


LH 


































C, 3LL 


all 


0 


0 


0 


0 


0 


0 


0 ; 


0 


0 


0 


0 


0 


0 


0 


0 



[0163] Layer 0 contains all data not quantized away by quantizer 0. This would be luminance data only: all of 3LL; 
all but 4 bitplanes of 2HL, 2LH, 3HL, 3LH and 3HH; all but 5 bitplanes of 2HH and all but 6 bitplanes of 1 HL and 1 LH. 
Layer 1 contains all data not in layer 0 and not quantized away by quantizer 1. This would be luminance bitplanes 5 
for 1 HL and 1 LH, bitplane 4 for 2 HH, bitplane 3 for 3HL and 3LH; all 3LL chrominance; all but 5 bitplanes for chromi- 
nance 3HL and 31 H; and all but 6 bitplanes for chrominance 2HL, 2LH and 3HH. Finally, layer 15 would contain the 
LSB of 1LH chrominance. 

Rate Control with Multiple Layers and Tile-Parts 

[0164] There several well known techniques for rate control in compression systems. The simplist method is to pick 
a distortion for every "unit" compressed (a unit may be an 8x8 block in JPEG, a frame in a motion sequence, a tile of 
a single image, a subband of a tile in a wavelet coded image, etc.). If the distortion selected leads to a bitrate higher 
than the desired average bitrate, the distortion allowed is increased for new units as they are compressed. If the dis- 
tortion selected leads to a bit rate lower than the desired average bitrate, the distortion allowed is decreased for new 
units as they are compressed. 

[0165] A more complex method buffers the compressed data from some number of "units." The bitrate and/or dis- 
tortion for each unit at each distortion level is stored. Then the distortion to allow across all the units in the buffer is 
determined when the buffer is full. If the buffer is sufficient to contain the entire image, extremely high quality results 
can be obtained. In JPEG 2000, layers are designed to contain increments to quality. Thus, selecting a distortion can 
mean selecting the number of layers to use for each code block or tile. A complete description of this type of rate control 
Is in, David Taubman, "High Performance Scalable Image Compression with EBCOT," IEEE Transactions on Image 
Processing. 

[0166] There are several disadvantages to this process. One disadvantage is that a buffer memory for the entire 
codestream is required. A second disadvantage is that the latency (time until any of the codestream is output) is high. 
A third disadvantage is that the second pass could take large amount of time. 

[0167] To mitigate these problems, each tile of a JPEG 2000 codestream is encoded as described above with at 
least two layers. At the completion of encoding each tile, a number of packets (e.g., layer, resolution, precinct, tile- 
component) are output to the codestream as a complete tile-part. The remaining layers are stored in the buffer. A 
second pass through the remaining coded data in the buffer is optional. During this second pass, extra packets from 
each tile are appended to the codestream as complete tile-parts as space or time allows. If in a fixed-rate application, 
then only packets within the given rate are appended. If in a fixed time application, then only number of cycles allowed. 
One embodiment of this process is shown in Figure 15A. Thus, these can be the 2 complete tile-parts outputfor each tile. 
[0168] Figure 15B illustrates a number of layers, layers 1-n. Layer 1 is output on the first pass, and the remaining 
layers are most likely below fixed-time or fixed-rate time limits. Layer 2 may be output on a second pass within fixed- 
time or fixed-rate requirements while achieving similar distortion over all the components. 

[0169] The above process is advantageous in that it allows the buffer to store a fraction of the coded data, the first 
data can be output (transmitted or stored) sooner, and the second pass through the data can be faster because there 
is less data to process. Also less memory is required lor buffering. 

[0170] The criterion for selecting which packets go into the first set of tile-parts can be similar to any other rate control 
algorithm. In one embodiment, the rate of packets can be less than the desired average bitrate for the whole image. 
For example, if a final compressed bitstream at 2.0 bpp is desired, the first pass could place 1 .5 bpp for every tile in 
the codestream, and buffer 1 bpp for every tile. 

[0171] The second pass can select from the remaining data the packets to place in the second tile part of^ach tile. 
Thus, to obtain a 2.0 bpp average encoding, some tiles that had high distortion after the first pass could receive all the 
remaining data saved for the tile, while other tile parts which had low distortion after the first pass might not have any 
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additional data transmitted. 

Rate Control for Compressed Codestream Data 

5 [0172] Some rate control techniques described herein include rate control performed on a compressed codestream 
based on a request implemented by selecting some number of layers to keep in the codestream. A parser may be 
used to produce a new codestream which shows the bitrate based on layers. This bitrate is equal to or less than the 
bitrate specified by the request. 

[0173J The parser may use a data structure referred to herein as a "packet structure." Note that this data structure 
10 may be used for other purposes such as, for example, the versatile pocket data structure described below. In one 
embodiment, the packet structure includes a packet start pointer and packet length. It also contains a tile number, a 
resolution, a component, layer, and a precinct the packet belongs to. Finally, it also consists of a selection flag. This 
flag, when set to a predetermined value (e.g., 1), indicates if the packet is selected in the array for writing out to a new 
codestream. 

is [0174] In one embodiment, packets are read in sequential order from a codestream based on the progression order 
information indicated by the COD marker. 

[01 75] The number of bytes is computed based on the bitrate desired by the request. The number of bytes belonging 
to layer 0 is added up to a total. Then this total of bytes is compared with the number of bytes desired. If the total is 
less than the number of bytes desired, one additional layer is added to the total. The process continues until the total 
20 is equal to or greater than the number of bytes desired or all packets have been added. 

[0176] During the process, those packets which have been added to the total, are marked as selected by the selection 
flag in the structure. 

[0177] If the total is equal to the number of bytes desired, the addition process is stopped. If the total exceeds the 
number of bytes desired, the packets in the last layer added are subtracted from the total. This is done to guarantee 
25 that the bitrate is below the bitrate desired. Consequently, during the subtraction step, packets which have been sub- 
tracted from the total are marked unselected. 

[0178] In one embodiment, the related markers such as SOT, COD, PLT are updated according to the request. 
Packets are written to the new codestream. The packet structure may be created using the following: 

30 



35 



40 



45 



SO 



55 



BNSDCCID: <EP T233624AiJ_> 



22 



EP 1 233 624 A1 



typedef struct JPACK_ { /* packet structure */ 

int start; /* packet starting point V 

int length; /* packet length V 

unsigned short t; /* tile number the packet belongs to *l 
unsigned short r; /* resolution the packet belongs to * / 
unsigned short c; /* component the packet belongs to */ 
unsigned short 1; /* layer the packet belongs to V 
unsigned short p; /* precinct the packet belongs to */ 
unsigned char select; /* selection flag */ 
} Pack_t; 



/* Store packets from tp->tile[i].Size[j] array to the packet structure array 
*//* Layer progression (LRCP) order */ 

if(progression_order == 0){ 
j = 0; 

for(i=0;i<number_of_tile;i-M-)( 
m = 0; 

for(l=0;l<layer;l++){ 
for(r=0;r<resolution+l,T++){ 
f or(c=0;c<component;c++) { 
f or(p=0 / -p<precinct[r],T3++) { 

tp->pk[j] -start = tp->tile[i] . pointer [m]; 
tp->pk[j].length = tp->tile[i]Size[m]; 
totaljength += tp->tile[i].Si2e[m]; 

tp->pk[j].t = i; 
tp->pk[j).r = r; 



tp->pk[j].l = I; 
tp->pk[j].c = c; 
tp->pk[j].p = p; 
m++; 

) 

} 

} 

> 

num_packet{i] = m; 

} 
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Versatile Packet Data Structure 

[01 79] The same packet data structure described above can be used to facilitate other parsing options, once packets 
are read into the structure. 

5 [0180] For resolution parsing, the packets which are to be excluded are marked unselected. For example, given a 
4 resolution codestream, and a request is to produce a 3-resolution codestream, a parser marks all packets which 
belong to resolution 4 unselected. Then the newly produced codestream contains only packets from resolution 1 up 
to resolution 3. 

[0181] Similarly, for component parsing, progression conversion parsing, quality parsing can be performed step by 
io step processing the packets in the structure. 

[0182] The packet data structure can handle complex requests. For example, a request which requires the parser 
to produce a codestream which has a 3-resolution, 2-layer, and 1 -component codestream. 

Clipping After Each Inverse Transform 

15 

[0183] As a result of quantization performed on wavelet coefficients, the final decoded pixels are often outside of the 
original range of allowed pixels from the specified bit depth. Typically, these pixels are clipped to the original range so 
that further image processing or display devices can use the original bit depth. 

[0184] For example, an eight bit image has pixel values between 0 and 255, inclusive. After lossy compression is 
20 used, the decoded image may contain values like -5 and 256. To provide an eight bit output, these values are clipped 
to 0 and 255 respectively. This clipping procedure always reduces pixel wise distortion because the original image did 
not contain pixels outside of the clipping bounds. This procedure is well known and recommend by the JPEG 2000 
standard. 

[0185] In addition to the bounds on the final output samples, there are bounds on the values coefficients can assume 
25 at the various stages of the wavelet transform. Just as quanti7ation can changp thp final decoded samples to lie outside 
the original bounds, quantization can change the partially transformed wavelet coefficients to lie outside their original 
bounds. If these coefficients are clipped to their original bounds, distortion will decrease. 

[0186] For example, after a horizontal (one dimensional) 5-3 reversible transform as specified by JPEG 2000 with 8 
bit input samples, the maximum value of the low pass coefficient is +191, and the minimum possible value is - 191 . 
The high pass coefficient must be between -255 and 255 inclusive. After the vertical one dimensional transform, the 
Low-Low coefficients are bounded by -286 and 287. Thus when decoding an eight bit image, when the first level low- 
low pass coefficients are generated (by the inverse wavelet transfrom from a higher level), the coefficients can be 
clipped to -286 and +287, and distortion will decrease. Likewise after the first level vertical inverse transformation is 
done, the low pass coefficients can be clipped to - 1 91 , +1 91 , and the high pass coefficients can be clipped to -255, 255. 

« [0187] For each subband, each filter, each transform level, and each image depth, there is a different maximum and 
minimum value for the coefficients. These maximum and minimum values can be computed by finding the signal that 
leads to the maximum and minimum and running the forward compression system and recording the maxima. The 
signals that lead to extreme values come from inputs where each pixel is either a maximum or minimum. Which pixels 
should be maximum and which pixels should be minimum can be determined by convolving sequences which are -1 

>o when the wavelet coefficient is negative and +1 when the wavelet coefficient is negative. Forthe 5-3 filter used in JPEG 
2000 Part I, the low pass signal of interest is [-1+1+1+1-1] and the high pass signal is [-1+1-1]. 
[0188] The signal (image) which will generate the largest LL value is: 
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(where +1 must be replaced by the input maximum (e.g., 255) and -1 must be replaced by the input minimum (e.g., 0). 
[01 89] For irreversible filters, it is not necessary to actually run the system to determine the maxima, simply convolving 
the wavelet coefficients is sufficient. For the reversible 5-3 filter, however the floor function is used in the computation 
of coefficients and is also used to determine the correct maxima. 
55 [0190] Note that this may be used for other filters (e.g., a 9-7 filter). 

[0191] Figure 28 is a flow diagram of one embodiment of a process for applying an inverse transform with clipping 
on partially transformed coefficients. The process is performed by processing logic, which may comprise hardware (e. 
g., circuitry, dedicated logic, etc.), software (such as that which runs on a general purpose computer system or a 
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dedicated machine), or a combination of both. 

[0192] Referring to Figure 2B, processing logic applies a first level inverse transform to coefficients {processing block 
2801 ). Thereafter, processing logic clips the partially transformed coefficients to a predetermined range (processing 
block 2802). Next, processing logic applies a first level inverse transform to the clipped coefficients (processing block 
5 2803) and clips the partially transformed coefficients to a predetermined range (processing block 2804), which is dif- 
ferent than the range in processing block 2802. Again, processing logic applies a first level inverse transform to clipped 
coefficients (processing block 2805) and clips the partially transformed coefficients to still anotherpredetermined range 
(processing block 2806). 

10 Simplified Colors pace Handling 

[0193] A typical decoding process including color management is shown in Figure 25. Referring to Figure 25, a file 
with a file format (e.g., afile format described in the JPEG 2000 standard) containing a restricted ICC profile is provided 
to a decoding device. Decompression block 2501 decompresses the file by taking the codestream portion of the file 

'5 and performing context modeling, entropy decoding, and applying an inverse wavelet transform, but does not perform 
color space operations. If the codestream indicates the RCT or ICT component transform should be used to decode 
the codestream, these will be performed by block 2502. That is, inverse RCT/ICT block 2502 takes the components 
and the "RCT Y/N" indication (RCT if yes, ICT is no) and performs the specified inverse transform and provides (non- 
display) RGB pixels. (If specified by the syntax, inverse level shifting is also performed.) 

20 [0194] Finally, the ICC color profile from the file format along with information about the display device will be used 
to produce the output pixels. 

[0195] Inverse ICC block 2503 receives the (non-display) RGB pixels and the ICC profile and applies an inverse 
color space transform to provide display RGB pixels. 

[0196] Figure 26 illustrates one embodiment of a non-preferred camera encoder. Referring to Figure 26, a camera 
25 generates YCrCb pixels. A converter 2602 converts the YCrCb pixels to RGB pixels and provides those two a typical 
JPEG 2000 encoder. The encoder comprises a RCT to ICT converter 2603 followed by a compressor 2604. The com- 
pressor generates an ICC A for codestream. 

[0197] Figure 27 illustrates one embodiment of a simpler camera encoder. That is, instead of including RCT/ICT 
converter 2603 and compressor 2604, a simple camera encoder includes only compressor block 2702. Referring to 
30 Figure 27, a camera 2701 generates YCrCb pixels and provides them to compressor 2702. Compressor comprises a 
JPEG 2000 encoder without an RCT conversion and generates an ICC B codestream with RCT equaling 1 (with syntax 
signaling that the inverse RCT should be used on decoding). The relationship between ICC B and ICC A is given by the 
following equation: 

ICC B = ICC A ° YCrCb 1 ° RCT 

where 0 represents function composition. 

[0198] Restricted ICC profiles are "syntaxes" for functions on pixels. A camera will typically write the same profile 
40 for all images, so ICC B is computed offline, and copied into each output file. In a prior art system there must be HW 
for YCrCb* 1 and RCT/ICT which operates on every pixel. 

Coding 4:2:2 and 4:1 :1 Data as 4:4:4 Data with Quantization 

45 [0199] The JPEG 2000 standard is typically used to handling data in a 4:4:4 format. It is not capable of describing 
how to reconstruct data in 4:1:1 or 4:2:2 formats in a 4:4:4 format for output. In one embodiment, when encoding 4:1 : 
1 data, the encoder treats 1 HL, 1 LH and 1 HH coefficients as zero. When encoding 4:2:2 data, the encoder treats 1 
HL and 1 HH coefficients as zero. Thus, with all information in the extra subbands quantized to zero, a decoder is able 
to receive the codestream in a way it expects. In other words, the encoded data resembles 4:4:4 data that has been 

50 heavily quantized. 

File Order for Thumbnail, Monitor, Printer, and Full Resolution and Quality 

[0200] Multiple images at multiple resolutions are important in many image processing situations. Depending on the 
55 application, a user may want to select different images of different resolutions. For example, thumbnail images may 
be used as an index into a large number of images. Also, a screen resolution image may be the image used to send 
to a monitor for display thereon. A print resolution image may be of lower quality for printer applications. 
[0201] In one embodiment, a codestream of an image is organized into sections so that different versions of the 
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image, such as, for example, a thumbnail version, a screen version, a print version and a lossless version, isprogressive 
by quality. 

[0202] In one embodiment, the packets are arranged such that certain packets correspond to particular resolutions 
such as a thumbnail. The combination of these packets with other packets represents the monitor resolution image, 

5 which when combined with other packets may represent the printer version, etc. Using the POC and tile parts, portions 
of a codestream may be grouped together. For example, all the tiles of the thumbnail size may be grouped together 
followed by tiles for another resolution followed by tiles of another resolution, etc. Figure 21 illustrates an example 
progression with tile parts for a single server. Each tile's thumbnail is grouped in tile-parts at the beginning of a file. 
Figure 21 A illustrates that tile-part 2101 is the only portion that is used for a thumbnail image. Figure 21 B illustrates 

10 that for a monitor resolution, tile-parts 2102-2104 have been included with tile-part 2101 . Figure 21C illustrates that 
for a printer resolution, tile-parts 2105 and 2106 have been included with tile-parts 2101-2104. Lastly, Figure 21D 
illustrates that for a lossless version of the data, the remaining three tile-parts 2107-2108 are included with the rest of 
the tile-parts. These sets of tile parts may be placed on a server in this progressive order. 

[0203] One embodiment of the process for accessing the groupings of tile parts is shown in Figure 1 6. The process 
*5 may be performed by processing logic that may comprise hardware (e.g., dedicated logic, circuitry, etc.), software 
(such as is run on a general purpose computer system or a dedicated machine), or a combination of both. The following 
steps assume that the image has been transformed with sufficient resolution levels and layers to divide the image into 
the four sizes. 

[0204] Referring to Figure 16, processing logic initially determines the correct resolution and layering for the thumbnail 
20 (processing block 1601). In one embodiment, to determine the correct resolution and layering for the thumbnail, 
processing logic creates a POC constrained to that resolution and layer for each tile and then creates a set of tile-parts 
and places this POC for each tile in the codestream. 

[0205] Next, processing logic repeats processing block 1601 for the monitor resolution given that the thumbnail 
packets arc already in the codestream (processing block 1602). Then, processing logic repeats processing block 1601 
25 for the printer resolution given that the monitor packets are already in the codestream (processing block 1603). 

[0206] Lastly, processing logic creates a POC marker with the extremes of the resolutions and layers for each tile 
(processing block 1604). In one embodiment, creating the POC with the extremes of the resolutions and layers is 
performed by creating a fourth set of tile-parts with the remaining tile-parts for a lossless version. 
[0207] Note that the particular orders of the packets defined in the POCs are not of importance, only the limits. 
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An Exemplary Computer System 



[0208] Figure 20 is a block diagram of an exemplary computer system. Referring to Figure 20, computer system 
2000 may comprise an exemplary client 150 or server 100 computer system. Computer system 2000 comprises a 
35 communication mechanism or bus 2011 for communicating information, and a processor 2012 coupled with bus 2011 
for processing information. Processor 2012 includes a microprocessor, but is not limited to a microprocessor, such as, 
for example, Pentium™, PowerPC™, Alpha™, etc. 

[0209] System 2000 further comprises a random access memory (RAM), or other dynamic storage device 2004 
(referred to as main memory) coupled to bus 2011 for storing information and instructions to be executed by processor 
40 2012. Main memory 2004 also may be used for storing temporary variables or other intermediate information during 
execution of instructions by processor 2012. 

[0210] Computer system 2000 also comprises a read only memory (ROM) and/or other static storage device 2006 
coupled to bus 2011 tor storing static information and instructions for processor 201 2, and a data storage device 2007, 
such as a magnetic disk or optical disk and its corresponding disk drive. Data storage device 2007 is coupled to bus 

45 2011 for storing information and instructions. 

[0211] Computer system 2000 may further be coupled to a display device 2021 , such as a cathode ray tube (CRT) 
or liquid crystal display (LCD), coupled to bus 2011 for displaying information to a computer user. An alphanumeric 
input device 2022, including alphanumeric and other keys, may also be coupled to bus 2011 for communicating infor- 
mation and command selections to processor 2012. An additional user input device is cursor control 2023, such as a 

so mouse, trackball, trackpad, stylus, or cursor direction keys : coupled to bus 2011 for communicating direction information 
and command selections to processor 2012, and for controlling cursor movement on display 2021 . 
[0212] Another device that may be coupled to bus 2011 is hard copy device 2024 : which may be used for printing 
instructions, data, or other information on a medium such as paper, film, or similar types of media. Furthermore, a 
sound recording and playback device, such as a speaker and/or microphone may optionally be coupled to bus 2011 

55 for audio interlacing with computer system 2000. Another device that may be coupled to bus 2011 is a wired/wireless 
communication capability 2025 to communication to a phone or handheld palm device. 

[0213] Note that any or all of the components of system 2000 and associated hardware may be used in the present 
invention. However, it can be appreciated that other configurations of the computer system may include some or ail of 
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the devices. 

[0214] Whereas many alterations and modifications of the present invention will no doubt become apparent to a 
person of ordinary skill in the art after having read the foregoing description, it is to be understood that any particular 
embodiment shown and described by way of illustration is in no way intended to be considered limiting. Therefore, 
references to details of various embodiments are not intended to limit the scope of the claims which in themselves 
recite only those features regarded as essential to the invention. 
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Claims 

1. A system comprising: 

a memory sized to include lines to store a band of an image and additional lines; 
a wavelet processing logic comprising 



a wavelet transform to generate coefficients when applied to data in the memory; 

access logic to read data from the memory into the line buffers to supply data stored in the memory to the 
wavelet transform and to store coefficients in the memory, such that after data stored at a first pair of lines 
is read from memory into the buffers of the access logic, the access logic reuses the first pair of lines to 
20 store coefficients generated by the wavelet transform that are associated with a second pair of lines dif- 

ferent from the first pair of lines. 

2. The system defined in Claim 1 wherein the access logic stores coefficients in contiguous lines of memory with 
coefficients from the same subband and decomposition level adjacent each other. 

25 

3. The system defined in Claim 1 wherein a first line of each of the first and second pairs of lines are located in the 
memory at an offset with respect to each other. 

4. The system defined in Claim 3 wherein the access logic stores the first outputs of the wavelet transform for each 
30 coefficient level in the additional lines within a distance of the offset. 

5. The system defined in Claim 3 wherein size of the offset is different for each transform level. 

6. The system defined in Claim 3 wherein the size of the offset is equal to: 

35 2< trans,orrn tevel of coefficient being stored) 

7. The system defined in Claim 6 wherein, during decomposition, the offset for storing the first rows of each pair of 
rows of L1 coefficients in the memory is two lines from the first row of data of the image associated with said each 
pair of rows of the L1 coefficients, and the offset for storing the first row of each pair of rows of L2 coefficients if 

40 four lines from the first row of L1 coefficients associated with said each pair of rows of the L2 coefficients. 

8. The system defined in Claim 1 wherein the access logic stores coefficients associated with a decomposition level 
greater than level three in the lines of the memory that previously stored the band of the image. 

9. The system defined in Claim 3 wherein the addition lines relating to the offset are above the line storing the band 
of the image. 

10. The system defined in Claim 1 wherein the wavelet transform is a forward wavelet transform. 
so 11. The system defined in Claim 1 wherein the wavelet transform is an inverse wavelet transform. 

12. A method comprising: 

reading data from a memory into line buffers to apply a wavelet transform thereto; and 
55 storing coefficients created by applying the wavelet transform at lines in the memory so that each set of -co- 

efficients generated from data stored at each pair of lines in the memory is stored in the memory at lines that 
are at an offset with request to said each pair of lines in the memory. 
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13. The method defined in Claim 12 further comprising access logic reusing a first pair of lines to store coefficients 
generated by a wavelet transform, that are associated with a second pair of lines different from the first pair of 
lines, after data stored at a first pair of lines is read from memory into the buffers of the access logic, and wherein 
a first line of each of the first and second pairs of lines are located in the memory at an offset with respect to each 

5 other. 

1 4. The method defined in Claim 1 3 further comprising the access logic storing the first outputs of the wavelet transform 
for each coefficient level in additional lines within a distance of the offset. 

io 15. The method defined in Claim 13 wherein size of the offset is different for each transform level. 

16. The method defined in Claim 13 wherein the size of the offset is equal to: 

2(transtorm level of coefficient being stored) 

is 17. The method defined in Claim 16 wherein, during decomposition, the offset for storing the first rows of each pair of 
rows of Li coefficients in the memory is two lines from the first row of data of the image associated with said each 
pair of rows of the L1 coefficients, and the offset for storing the first row of each pair of rows of L2 coefficients if 
four lines from the first row of L1 coefficients associated with said each pair of rows of the L2 coefficients. 

20 18. The method defined in Claim 12 further comprising access logic storing coefficients associated with a decompo- 
sition level greater than level three in the lines of the memory that previously stored the band of the image. 

19. The method defined in Claim 13 wherein the addition lines relating to the offset are above the line storing the band 
of the image. 

25 

20. An article of manufacture comprising at least one recordable media storing executable instructions thereon which, 
when executed by a processing device, cause the processing device to: 

read data from a memory into line buffers to appfy a wavelet transform thereto; and 
30 store coefficients created by applying the wavelet transform at lines in the memory so that each set of coeffi- 

cients generated from data stored at each pair of lines in the memory is stored in the memory at lines that are 
at an offset with request to said each pair of lines in the memory. 

21 . The article of manufacture defined in Claim 20 further comprising instructions, which when executed by the process- 
es ing device cause the processing device to reuse a first pair of lines to store coefficients generated by a wavelet 

transform, that are associated with a second pair of lines different from the first pair of lines, after data stored at a 
first pair of lines is read from memory into the buffers of the access logic, and wherein a first line of each of the 
first and second pairs of lines are located in the memory at an offset with respect to each other. 

^0 22. An apparatus comprising: 

means for reading data from a memory into line buffers to apply a wavelet transform thereto; and 
means for storing coefficients created by applying the wavelet transform at lines in the memory so that each 
set of coefficients generated from data stored at each pair of lines in the memory is stored in the memory at 
45 lines that are at an offset with request to said each pair of lines in the memory. 
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FIG.1 8 
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